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Abstract 



Two new QCD sum rules for nucleon tensor charge are derived from a 
mixed correlator of spin-1/2 and spin-3/2 nucleon interpolating fields. These 
sum rules are analyzed along with a sum rule obtained from the usual cor- 
relator of a general spin-1/2 nucleon interpolating field. The validity and 
reliability of the sum rules are examined by monitoring the contaminations 
arising from the transitions and continuum and the convergence of the oper- 
ator product expansion. Valid sum rules are identified and their predictions 
are presented. It is found that the vacuum tensor susceptibility induced by 
the external field plays an important role in determining both the validity and 
predictions of the sum rules. The uncertainties associated with the sum-rule 
predictions are also discussed. 
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I. INTRODUCTION 



Tensor charge of the nucleon, defined through the nucleon matrix element of the tensor 
current, qaapQ, at zero momentum transfer, is perhaps the least known nucleon charge both 
theoretically and experimentally. Like other charges of the nucleon, the tensor charge reveals 
important information on the nonperturbative structure of the nucleon and understanding 
it from the underlying theory of strong interaction, quantum chromodynamics (QCD), is of 
great importance. 

The nucleon tensor charge 6q is defined by 

{N\qa^pq\N) = 6qUa^pU , (1) 

where U stands for the nucleon Dirac spinor. It has been shown by Jaffe and Ji [jl| that the 
first moment of the twist-two transversity quark distribution in a nucleon, hi{x), is related 
to the nucleon tensor charge 







1 _ 

dx h\{x) — h^{x) = 6q , (2) 



where hi and receive contributions from quark and antiquark (of fiavor g), respectively. 
Unlike the twist-two spin averaged quark distribution fi{x) and helicity difference quark 
distribution gi{x), the hi{x) flips chirality and is thus suppressed in inclusive deep inelastic 
lepton scattering p|. As such, there is no existing experimental data on hi{x) and the 
nucleon tensor charge. However, there is no such suppression of hi{x) in deep inelastic 
processes with hadronic initial states such as Drell-Yan [^0, and experiments have been 
planned to measure the quark transversity distribution in the nucleon at RHIC [Q], HERA 
§, and CERN g. 

Recently, the nucleon tensor charge has been investigated theoretically in various models 
and approaches |]7|-|TI|. Previous studies of the nucleon tensor charge via the QCD sum- 



rule method have been made by He and Ji These authors have found that the sum 



rules obtained from the correlator of loffe's nucleon interpolating fleld |T2| are unstable, 
making the extraction of the nucleon tensor charge from those sum rules difficult. In this 
paper, we derive two new QCD sum rules from the mixed correlator of spin-1/2 and spin- 
3/2 interpolating flelds using the external-fleld method. The consideration of the mixed 
correlator is motivated by its success in nucleon mass sum rules |]TB[ and in the sum rules 



for nucleon axial- vector coupling constants |14]. We analyze these new sum rules along with 



a sum rule obtained from the usual correlator with a general spin-1/2 nucleon interpolating 
field. We examine the performances and reveal the validity and reliability of these sum rules. 

A well-known problem that arises in calculating hadronic matrix elements within the 
QCD sum-rule method is the unwanted physics associated with transitions from ground state 
to excited states. The contributions of these transitions are not exponentially suppressed 
relative to the ground state contribution which contains the ground state property of interest 



T5| . In practice, there are two formalisms for treating the transitions. The usual formalism 
is to model the sum over all the transitions from ground state to excited states by a single 
unknown parameter. This parameter is extracted from the sum rules, along with the ground 
state property, and the continuum threshold is usually fixed at the same value as that for 
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the corresponding mass sum rule. The resolution of the ground state signal relies on the 
different polynomial (in Borel mass) behaviors of the ground state signal and the transition 
term. 

An alternative formalism for treating the transitions is to multiply the invariant functions 
by the factor M| — (with Mh the ground state hadron mass) before performing the Borel 
transform. The transition contribution is then exponentially suppressed relative to the 
ground state contribution and can be included in the continuum model in practice. In 
this case, it is important to treat the continuum threshold as an unknown parameter to be 
extracted from the sum rule along with the ground state property. Assuming the threshold 
to be the same as that for the mass sum rule is unjustified and may introduce artificial bias 
to the extracted hadron properties fl^. We observe that in the previous works 0,11 the 
alternative formalism was adopted but the continuum threshold was taken to be the same 
as the one used in the nucleon mass sum rule. In the present paper, we shall consider both 
formalisms. 

Both sides of QCD sum rules depend on an auxiliary parameter-Borel mass, which is 
introduced through the use of the Borel transform. If a sum rule were perfect, one would 
expect that the two sides of the sum rule overlap for all values of the Borel mass. However, 
in practice, the two sides of the sum rules overlap only in a limited range of the Borel mass 
(at best) because of the truncation of the operator product expansion (OPE) and the crudity 
of the models for the transitions and continuum. Thus, in order to extract the properties 
of the ground state by matching the sum rules, one should work in a region of Borel mass 
where the ground state contribution dominates the phenomenological side. This usually 
sets an upper bound in the Borel mass space, beyond which the excited-state and transition 
contributions dominate and the background noise is hence too strong to reliably isolate the 
ground state signal. On the other hand, the truncated OPE must be sufficiently convergent 
as to accurately describe the true OPE. This, in practice, sets a lower limit in the Borel 
mass space, beyond which higher order terms not present in the truncated OPE may be out 
of control. 

Therefore, we expect a sum rule to work if the two sides of the sum rule match in a valid 
window in Borel mass space where the ground state contribution of interest dominates the 
phenomenological side and the higher order OPE terms are under control. The validity and 
reliability of a sum rule is then determined by the quality of overlap of the two sides, the size 
of valid Borel window, and the relative contributions of the excited states and the highest 
OPE terms in the valid Borel window. A sum rule failing to have a valid Borel window 
should be considered invalid, and results obtained from such an invalid sum rule may be 
meaningless and misleading. In this paper, we will analyze the sum rules with respect to 
the above criteria. 

This paper is organized as follows. In Sec. ||, we derive QCD sum rules for the nucleon 
tensor charge within the external-field approach. In Sec. we analyze the sum rules and 
present results and discussions. Section |3 is devoted to summary and conclusion. 
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II. SUM RULES FOR NUCLEON TENSOR CHARGE 



In this section, we first give a brief outline of the external-field QCD sum-rule method 
with an emphasis on some issues related to the treatment of the transitions between ground 
state and excited states. We then use this method to derive QCD sum rules for the nucleon 
tensor charge. 



A. Outline of the method 



Consider the nucleon two-point correlator in the presence of a constant external (anti- 
symmetric) tensor field, Z"^, defined by 



n(g) =i / dW«-^(0|Tr7(x)r7(0)|0), 



(3) 



where 77 is a nucleon interpolating field, constructed from quark fields, which couples to the 
nucleon. To lowest dimension, there are two independent interpolating fields which couple to 
spin-1/2 states only and carry the quantum numbers of the nucleon. So, a general spin-1/2 
interpolating field for the proton can be expressed as 



where /3 is a real parameter and 

mix) = ^abc 



V2[X) 



^abc 



u''^{x)Cd\x)\ 75U'(X) . 



(4) 

(5) 
(6) 



Here u{x) and d{x) denote up and down quark operators, a, 6, and c are color indices, and C 
is the charge conjugation matrix. The interpolating field for the neutron can be obtained by 



changing u{d) into d{u). The interpolating field advocated by loffe [|T2| and used in previous 
sum-rule calculations of the nucleon tensor charge may be recovered by setting [3 = —1 
and multiplying an overall factor of —2. 

There is a spin-3/2 interpolating field which also couples to the proton through its spin- 
1/2 component pTIjH 



The couplings of 77^/^ and r]^^"^ to the proton are defined as 

{Q\V'/'\N) = \,U{q) , (iViTjf |0) = A2F(g)75 + 7. 



(7) 



(8) 



where U{q) denotes the Dirac spinor of the proton with the normalization U {q)U{q) = 2Mn, 
and Ai and A2 describe the coupling strengths of the interpolating fields to the proton. 

Therefore, independent QCD sum rules can be obtained from the correlators of various 
combinations of 7/^/^ and t/^^^. In the present work, we will consider the mixed correlator 
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of 7]^/^ and rfj^"^ as well as the 1/2 correlator of two r/^/^'s. While the latter, at /5 = — 1, 
leads to the sum rules discussed previously, the former gives rise to new independent sum 
rules for the nucleon tensor charge. In Ref |T^, it has been demonstrated that the nucleon 



mass sum rules from the mixed correlator have broader valid Borel windows and relatively 
smaller continuum contamination as compared to the sum rules from the 1/2 correlator. 

The subscript Z in Eq. (|^) denotes that n(g) is evaluated by adding a tensor coupling 
term 

^L = Z'^^Y.S,l{x)a^^q{x) , (9) 

to the QCD Lagrangian, where q = {u, d, s}. Here Qg keeps track of the flavor structure of 
the external field. For example, if one is interested in the isovector (isoscalar) tensor charge, 
g^ = 6u- 6d (5f^ = 6u + 6d), one simply sets gu = -Qd = I {Qu = Qd = 1) and Qs = 0. The 
external-field method proceeds by expanding Il{q) to first order in the external field 

n(g) = nW(g)+Z-^n2(g) + ---, (10) 



where Il^^\q) is the correlator in the absence of the external field from which the usual 



nucleon mass sum rules are derived. It is the linear response to the external field, Ill^l{q), 



that gives rise to the sum rules for the nucleon tensor charge. 

The QCD side of the sum rule is obtained by carrying out the OPE. The external field 
can interact with the quark field directly; it also polarizes the QCD vacuum and induces 
nonperturbative condensates. The latter are summarized by vacuum susceptibilities, which, 
to first order in the external field, are defined as 

{l(^apq)z = 9qX{l(l)ZaP , (H) 

{qgsGaf3q)z = 9qf^{Q(l)Zaf3 , (12) 
{q9sGaf3l5q)z = -Wq^Q(l)Zaf3 , (13) 

where gs is the strong coupling constant (Z)^ = + igsAp), Gap = \^ai3p\G^^ with Gap the 
gluon field tensor, x, k, and ^ are vacuum susceptibilities, and (O) denotes the usual vacuum 
condensate. Here, the isospin breaking effect has been neglected and will be henceforth. The 
strangeness contribution will also be ignored in the present paper as it is expected to be small 



PJTT[| . The calculation of the Wilson coefficients is straightforward following the techniques 
discussed extensively in the literature. Some details of the calculations are presented in the 
Appendix. 

We now turn to the phenomenological description, which can be obtained by expand- 
ing Ii^^^p{q) in terms of physical intermediate states. For a generic invariant function, the 
phenomenological representation can be written as 

_ ^o9t , Cn^n* 

^ {Ml - g2)2 + ^ (M^ - q^){Ml, - q^) 

-|-terms involving only excited states , (14) 
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where Mat and Mat* are the masses of the ground state nucleon and the excited state, 
respectively, and Xq = A1A2, for the mixed and 1/2 correlators, respectively. The first 
term contains the nucleon tensor charge g^j, = J2q9q^Q of interest, the second term is the 
transition term from ground state to excited states, and the rest involves only the excited 
states. After a usual Borel transformation, one finds 



Cn^N* _ (M^,_M^)/M= 



-Ml/M^ 



+exponentially suppressed terms , (15) 

where M is the Borel mass. It can be seen clearly that the second term is not exponentially 
damped as compared to the first term. The usual formalism is to approximate the second 
term as Ae~*^jv/^ with A a phenomenological constant to be extracted from the sum rule 
along with g^. 

An alternative formalism is to use the combination [Mj^ — q^] Y\.^^\q^)^ instead of n''^^(g^) 
alone. The phenomenological side then becomes 

n(i) (M^) ~ A^^^e-*^^/*^' + CN^N*e-^'^N'/^^ + exponentially suppressed terms . (16) 

N* 

The transitions from ground state to excited states (second term) are now exponentially 
suppressed with respect to the ground state (first) term and hence can be absorbed into 
the continuum model. As such, one no longer needs to introduce any phenomenological 
parameter to represent the transitions from the ground state to excited states as their 



contribution has been included in the continuum model. However, as stressed in Ref. |T6 
the continuum threshold must be treated as an unknown parameter to be extracted from 
the sum rule. 

Since the linear response n^^^j in general contains distinct Dirac and Lorentz structures, 
one may obtain many sum rules, one for each invariant structure. However, these sum 
rules do not work equally well in practice. In particular, some sum rules work well while 
the others may fail. Recently, it has been pointed out by the present authors that 
chirality plays an important role in determining the reliability of baryon sum rules. There it 
was argued that for light baryons the chiral-odd sum rules (where the chiral-odd operators 
dominate) are generally more reliable than the chiral-even sum rules (where the chiral-even 
operators dominates). Thus, we will only consider the chiral-odd sum rules and disregard 
the chiral-even sum rules in the discussions to follow. 



B. Sum rules from the mixed correlator 

The mixed correlator is defined by 



H£)(g)+Z-X^) (g) + -- - , (17) 
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where -qj/"^ = 7^75?]^''^ and r/^/^ and 77^/^ are given in Eqs. @ and (0), respectively. The use 
of 7]j/'^ is to facihtate the calculations and does not affect the sum rules. The linear response 
of n^iy(g) to the external field can be decomposed into various invariant structures: 



n 



+n3(g2) _ q^^^) qq^ + --- , 



(18) 



where q = 5^7'^. Here, the first three structures receive contributions from the ground state 
proton and lead to chiral-odd sum rules. Since there is no continuum model for the sum rule 
from lis, we only focus on the sum rules from IIi and 112. Adopting the usual formalism for 
the transitions, we obtain the sum rule from IIi 



+ 



24 

ab 



C4 



^8/27 



rriQa 
~2A 



K,a 



C2 



L-^'"' In ^ 



C5 



^8/9 




7e 



+ 



"^0^ „ r-2/27 



96 



C3^ 



864M2 
and the sum rule from 112 
2a 



144M2 ^ 72M2 



^8/9^ 



1728M2 



C7 



^-6/27 



Vm2 



(19) 



rflEoi^4/9^2_!!^^^-2/27 
3 1 32 



In 



'M2 



-7e 



^^^8/27 ^ ^ ^ . 8/9 ^ ^ r 8/9 



Xmga2 



9M2 



ah 



216M2 



= 4 



A1A2 



M, 



iV 



9M2 
M2 




^ , ^-2/27 



Mjr/M'^ 



(20) 



where a = -(2^)2(gg), h = {glG^), ml = -{qg,a ■ Gq)/{qq), \, = (27r)2Ai, A2 = (27r)2A2, 
Eq = 1 — e^^o/*"^ with So the continuum threshold, and 7^,^ = 0.577- ■ ■ the Euler-Mascheroni 
constant. The anomalous dimensions of various operators have been taken into account 
through the factor L = In^M"^ / Aq^^)) / 1ii(/^^/Aqcd); where /x = SOOMeV is the renormaliza- 
tion scale and Aqcd = ISOMeV is the QCD scale parameter. Here A denotes the infrared 
cut-off, which arises from the quark propagator in the presence of both the external field 
and background gluon field (see Appendix); a reasonable choice for A is A 
and dj are defined as 



19 . The Ci 
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ci = (5 + 9P)gu + (1 - 9P)gd , 

C3 = -(29 + 85P)gu + (5 - 7ip)gd , 

C5 = 69(3gu + (-6 + 2(3)gd , 

c, = -(47 + 87P)g^ + (23 - 237/3)^, 

di = {3 + AI3)gu + 5/3^rf , 

4= (2 + 13/3)5„ + (l + ll/3)^d , 

4 = (10-13/3)^„-(4 + 2/3)^, , 

^7 = -(158 + 10/3)^„ + (7 + 71/3)5<i , 



C2^{1- (5){^gu- Qd) , 
C4 = -7il + P)gu + il-np)gd 
C6 = (7 + 9/3)^„- (1 + 21/3)5,, 
cg = 5(1 + -(8- 28/3)5, , 

c?2 = (1 - /3)(fi'« - fi-d) , 

d4 = (7 + 3/3)5„- (1-3/3)5,, 

de = 2/35„ + 6/35, , 

4 = (5 + 7/3)5„ + (1 + 23/3)5, • 



(21) 



When the alternative formahsm for the transitions is used, the sum rule from Hi is 



96 

' 72 

ab 
864 



C6 1 + 



C8 1 + 



^8/9 ^ X^ga' . f-, , Mj, 



M2 y 



1728 



C7 1 + 



Ml 
M2 



jL^/9^M^,AiA2 5Te" 



^-6/27 



(22) 



and the sum rule from 112 is 
2a 



^-2/27 



6 



M 



N 



M2 



^8/9 



9 
a6 



1 + ^^ ^8/9_X^^^^^ 



M2 



216 



) r -6/27 
M2 i 



216 



4 1 + 



M 



N 



M2 



= 8MjvAiA2 5^6 



(23) 



where = 1 - e-"o/*^'(so/M2 + 1). 



C. Sum rules from the 1/2 correlator 

The 1/2 correlator is given by 
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n(g) = z / d^xe"'''{0\T[r]'/\x)t/\0)]\0), 



= nW(g) + Z"^n«(g) + -- - . 
The linear response term in this case has three distinct structures 



n 



(1) 

a/3 



n4(g^) (gcrQ/3 + o-Q/sg) + n5(g^) aa/s + ^eiq^) i {qalp - qpla) q ■ 



(24) 



(25) 



However, only 114 leads to a chiral-odd sum rule, which is given by 



32^1 ° 96 



In 



A2 



-7e 



2 

^-10/27 _ e ^-10/27 

192 



r I ^'^ /. r 16/27 I r T 16/27 

48 ■' 576M2 96M2 ■'^ 1728M2 



^^^-14/27 



a6 



6912M2 
in the usual formalism and 
a 



(26) 



< 1^ I A2 



^-10/27 



"T92 



+ 



48 
a6 



/6^^^/^^ 1 + 



48 










1728 


M2 i 





'^^^ . r 16/27 /^l , 

576 ^^^ 



M2 



(l + 



6912"" V 
in the alternative formalism, where 



(27) 



h = 


(10 + 6P^)g^ - (3 + 2/3 - /32)(?, , 


/2 


= (1 - /?) V , 




h = 


{A-10P + 6P^)gu-5il-P^)gd, 


/4 


= -2{l-f3)gu+{l 




h = 


(20 + 2/5 - 22/32)(7„ + 9(1 - P')gd , 


/6 


= 2{1-P)g„-{1- 




f7 = 


(6-7/5 + /32)(7„-16(l-/32)(7d, 








/8 = 


(62 + 36/3 + 70f3^)gu + (1 - 6/? + 17/5^)^^ . 









(28) 

The sum rule Eq. (P7|), at /5 = — 1, can be compared with Eq. (20) of Ref. [Q. While the 
first and fourth terms agree, the rest disagree. The factor ln(M2/A2) — M^7g^ — M'^Eq 
in the second term has the form M^ln(M^/A2) — in Ref. [||. We have also compared 
the corresponding form of this term before the Borel transformation and found an overall 
sign difference. Furthermore, the Borel transformation of this term has not been done 
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correctly in Ref. |^ (see Appendix for more details). The third term is absent in Ref. [§. 
The differences in the last four terms arise from the incomplete calculation of the Wilson 
coefficients in Ref. ^ (see Appendix). There is also a minor difference in the anomalous 
dimension for the mixed quark-gluon operator qgga-Gq. We use —2/27, while Ref. |^ simply 
takes 0. 



III. SUM-RULE ANALYSIS 

We now analyze the sum rules derived in the previous section. To this end, we optimize 
the fit of the two sides of each individual sum rule in a valid Borel window, which is chosen 
such that the highest order (in 1/M^) term contributes no more than ~ 10% to the QCD 
side while the contribution of the transition plus continuum is less than 50% of the total 
phenomenological side (i.e., the sum of the contributions from ground state, transition, and 
continuum). Those sum rules which fail to have a valid Borel window is considered invalid 



and useless. Such criteria have been adopted in various sum-rule calculations [ll3| , p!^ , |20 
While the selection of 50% is obvious for the ground state dominance, the selection of 10% 
is a reasonably conservative criterion. To quantify the fit of the two sides, we sample the sum 
rule in the valid Borel window and minimize the measure (5(M^) = (LHS-RHS)^ averaged 
over 50 points evenly spaced in the valid Borel window. Here, LHS and RHS denote the 
left- and right-hand sides of the sum rule, respectively. 



A. Inputs 

The inputs for the sum rules can be classified into three categories: a) the usual vacuum 
condensates; b) the vacuum susceptibilities and c) the nucleon mass Mat, the coupling 
strengths Ai and A2, and the continuum threshold sq. [The continuum threshold is needed 
as input for the sum rules (^^-^^) and (^B|).] For the values of vacuum condensates we use 
the central values given in Ref. a = 0.52 GeV^, h = 1.2 GeV^ and = 0.72 GeV^. 



The mixing parameter P can be chosen to minimize the overlap with the excited states and 
broaden the valid Borel window. Here we adopt (3 = for the sum rules from the mixed 
correlator and (3 = —1.2 for the sum rules from the 1/2 correlator, as determined in Ref. |]13 . 

The vacuum susceptibilities arise from the response of condensates to the external field. 
It is easy to show that to first order in the external field 

(qOcfsq), = Z'<' Wm i [ d^xe'^'^ {T[Y^ g,q{y)a,sq{y), ^./^(O)]) , (29) 

where C^/s = {qo'aisq, qgsGafsq, qdsG'j^q}. The two-point correlation function in Eq. (|29D 
can be decomposed into 

i / d^xe'Py{T[^ggq{y)a^sqiy),Oa/3iO)]) = {g^^gsp - g^p g5c.)^i{p^) 

+ {g-ic PSPP - 95a - 5-7/3 PsPa + 95/3 PjPa) ^2 (P^) • (30) 
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Neglecting isospin breaking effect and strangeness contribution and using Eqs. (pl]-|T3|), we 
find 



Xa = -Stc^U^O) , Ka = -Stt'^UKO) , = -SnHufiO) . (31) 

These definitions are equivalent to those given in Ref. P). It is worthwhile emphasizing 
that the perturbative contribution to (0) has to be subtracted as it has already been 
included through the perturbative quark propagator in the presence of the external field 
(see Appendix). 

Obviously, it is difficult to evaluate Hf (0) directly from QCD at this stage. The QCD 
sum-rule method has been invoked to estimate Hf (0) previously. The results given in Ref. |§] 
are ~ — O.lSGeV^ and = ~ O.lGeV^. Recently, the authors of Ref. have argued 
that the estimate of the vacuum tensor susceptibility x R-^f. is very rough, and their 
detailed analyses lead to a value for x which has opposite sign and much larger (~ 4 times 
larger) magnitude. Given this uncertainty, we will consider the x values in the following 
range 

- 0.5 GeV^ < xa < 1-5 GeV^ . (32) 

Since k and ^ are small and the sum rules are relatively insensitive to their values, we use 
the values of k and ^ given in Ref. |^ for simplicity. 

The extraction of the nucleon tensor charge from the sum rules also requires the knowl- 
edge of the nucleon mass Mj^, the couplings Ai and A2, and the continuum threshold Sq 
(for the sum rules in the usual formalism). These parameters can be determined from the 
corresponding mass sum rules. Here we take the central values obtained in Ref. [|T^| 



Mn = 0.96GeV , A1A2 = 0.41GeV'' , sq = 1.69GeV' , (33) 
for the sum rules from the mixed correlator and 

Mn = l.lTGeV , A^ = 0.20GeV^ , sq = 2.34GeV2 , (34) 



for the sum rules from the 1/2 correlator. Note that the nucleon mass in Eq. (|3^) is very 
close to its experimental value. It is, however, much larger than the experimental value in 
Eq. (|3^). In the literature it is often found that the nucleon mass is fixed at its experimental 
value and the mass sum rules are used to determine the couplings and continuum threshold. 
This may introduce artificial errors. Moreover, when the same criteria are applied to the 
analyses of the nucleon mass sum rules, we find that the continuum threshold cannot be 
searched successfully when the nucleon mass is held fixed at its experimental value. 



B. Results 



Let us start from the sum rules in the usual formalism, Eqs.(TP-|5D|) and (^). The 
nucleon tensor charge (7^ and the transition strength are extracted from the sum rules, 
with the continuum threshold fixed at its value for the corresponding nucleon mass sum 
rule. We find that for the x values interested here there does not exist valid Borel window 
for the sum rule Eq. (EDI) according to the criteria set above; it is thus invalid. 
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FIG. 1. Sum-rule predictions for the isovector nucleon tensor charge gj^ as functions of x- 
The sohd and dashed curves correspond to the sum rules Eqs. ( [l9| ) and (|26|), respectively. The 
dotted curve is also from the sum rule Eq. (p6|) but with = 0.94GeV, Af = 0.26GeV^, and 
So = 2.3GeV2 as given in Ref. §. 

The sum rule Eq. (|^) works only for a limited range of x values. The predicted isovector 
((7^) and isoscalar (g^) nucleon tensor charges as functions of x are plotted in Figs. |I| and 
0, respectively. Both and g^ increase as x increases with essentially the same rate. This 
indicates that 6u strongly depends on x while 5d has much weaker dependence on x- One 
can see that g^ > g^ > 0, implying 6u > 6d > 0. We also find that the size of the valid 
Borel window shrinks with decreasing x- 

The sum rule Eq. (^) is valid for all the x values considered here. Its predictions for the 
nucleon tensor charges are also given in Figs. and ^. The predicted tensor charges increase 
as X increases with a similar rate as those obtained from Eq. (JT^. The predictions again 
show 6u > Sd > 0, but 6d has stronger x dependence than that found from Eq. (|l^). For a 
given X value, the size of the valid Borel window is in general larger than that for Eq. ([T^. 

One notices that the predictions from Eq. (^) are significantly larger than those from 
Eq. (|TUp. This partially attributes to the large nucleon mass and small coupling [see Eq. 
used in Eq. (p6|). To further illustrate this point, we also displayed in Figs. ^ and |^ the 
predictions obtained from the sum rule Eq. (^), with = 0.94GeV, = 0.26GeV^, and 
So = 2.3GeV^ as quoted in Ref. 0. In this case, the sum rule Eq. (pB]) only works for a 
limited range of x values and the results are much closer to those from Eq. ([191) . 

We now turn to the sum rules in the alternative formalism, Eqs. (p2|-p^) and (pT]). As 
emphasized earlier, the continuum threshold Sq must be treated as an unknown parameter 
to be extracted from the sum rules along with the tensor charges. We find that the sum 
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FIG. 2. Sum-rule predictions for the isoscalar nucleon tensor charge as functions of x- 
The sohd and dashed curves correspond to the sum rules Eqs. ( [l^ ) and (|26|), respectively. The 
dotted curve is also from the sum rule Eq. (p6|) but with = 0.94GeV, Af = 0.26GeV^, and 
So = 2.3GeV2 as given in Ref. §. 

rule Eq. ([23| ) fails to have a valid Borel window and is thus invalid. The sum rules Eqs. ( p2D 
and ( P7D are also invalid for the case of isoscalar tensor charge. For isovector tensor charge, 
the sum rule Eq. (|19]) only works for small x values (x < 0.1) and its prediction is negative 

(~ —1.02 0.16). These results are much smaller than the corresponding ones in the 

usual formalism. On the other hand, the sum rule Eq. (^) is only valid for large x value 
(x ^ 0.8) and the predictions are slightly larger than those found from Eq. (^61) . 

C. Discussion 

To see how well the valid sum rules work, we have plotted in Fig. ^ the left- and right- 
hand sides of the sum rules Eqs. (|T^ and (^) for x = 1-0, with and Ai extracted from 
the sum rules. The corresponding valid Borel windows and the relative contributions of the 
transition plus continuum and the highest order term in the OPE are displayed in Fig. ^. 
One notices that the sum rule Eq. ( p6D is valid in a broad Borel regime where the transition 
plus continuum contributes about 20% at the lower bound and the required 50% at the 
upper bound. The sum rule Eq. (0) has a much smaller valid Borel regime where the 
contributions from the transition and continuum are greater than 40% throughout the valid 
Borel window. On the other hand, the overlap of the two sides of Eq. (|T^) is better than 
that of Eq. (|26|). This pattern is seen at other x values. It is therefore hard to compare the 
reliabilities of these two sum rules directly. 
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FIG. 3. The left- and right-hand sides of the sum rules Eqs. (|T9|) and (^) as functions of the 
Borel mass for x = 1-0, with and Ai extracted from the sum rules. The four curves correspond 
to the left- (solid curve) and right-hand (dashed curve) sides of Eq. ( p^ ) and the left- (dotted curve) 
and right-hand (dot-dashed curve) sides of Eq. (p6|), respectively 



As shown in Eq. (|T5|), the contribution from the transitions is in general a complicated 
function of the Borel mass, which is crudely modeled by a constant {Ai) in the usual formal- 
ism. This will change the curvature of the phenomenological side as a function of the Borel 
mass and may lead to errors to the extracted nucleon tensor charges. It is unclear how to 
reliably estimate the size of these errors within the criteria discussed in the present paper, 
though some estimates have been given in the study of other nucleon properties , 22 . 

Moreover, the spectral parameters (mass, coupling, and continuum threshold) appearing 
in the nucleon mass sum rules also enter the sum rules for the nucleon tensor charge. This 
means that the uncertainties associated with these parameters will give rise to additional 
uncertainties, in addition to the uncertainties in the sum rules for the tensor charge them- 
selves. This is a general drawback of the external field sum-rule approach, and we have 
demonstrated this point explicitly here. Clearly, one needs accurate knowledge of the vac- 
uum spectral parameters in order to extract the tensor charges cleanly, which is difficult 
to achieve given the current implementation of the QCD sum-rule approach. Therefore, we 
expect large uncertainties associated with the predictions of the valid sum rules in the usual 
formalism presented above. 

In the alternative formalism, the contribution of the transitions is absorbed into the 
continuum. It has been argued recently by one of us that this alternative formalism 
may have the potential to improve the reliability of external-field sum-rule calculations, as 
the transitions are explicitly modeled by the modified continuum model. However, we have 
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found here that the sum rules in the alternative formalism fail in most cases interested 
here. This may indicate that the information contained in the truncated OPE is not enough 
and/or the transition and continuum give rise to too much background noise to isolate 
the information about the nucleon tensor charge. This might also suggest that the reason 
for the sum rules in the usual formalism to work is because the transition contribution is 
approximated by a constant. 

We have emphasized throughout this paper that in the alternative formalism the con- 
tinuum threshold has to be treated as a unknown parameter because the continuum model 
is modified relative to that in the usual formalism and contains the transition contribution. 
Fixing the continuum threshold at its value for the corresponding mass sum rule may in- 
troduce arbitrary artificial effect to the extracted tensor charges. In Ref. f^, the sum rule 
Eq. ( p7D is used in extracting the nucleon tensor charges but with the continuum threshold 
fixed at its value for the nucleon mass sum rule. However, we have found that this sum rule 
is invalid for xa = — O.lSGeV^ as used in Ref. @, and hence the results obtained from this 
sum rule are unreliable. 

In their analyses, the authors of Ref. |^ express the tensor charge as a function of the 
Borel mass by dividing both sides of Eq. (pT]) by MatA^ e~^^'v/*^ and find the tensor charge 
changes quickly with the Borel mass. This is actually an indication of the poor overlap 
of the two sides of the sum rule. To illustrate this point, we have plotted in Fig. ^ the 
left- and right-hand sides of Eq. with the use of Mjy = 0.94GeV, = 0.26GeV^ 

Xa = — O.lSGeV^, and the continuum threshold fixed at sq = 2.3GeV^ as quoted in Ref. [§. 




(GeV^) 



FIG. 4. Relative contributions of the transition plus continuum and the highest order OPE 
terms to the sum rules as functions of the Borel mass. The dashed and dot-dashed curves correspond 
to the sum rules Eqs. (|l^) and (p6|), respectively. 
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FIG. 5. The left- (solid curve) and right-hand (dashed curve) sides of the sum rule (^) as 
functions of the Borel mass, with the use of Mn = 0.94GeV, Af = 0.26GeV6, x« = -O.lSGeV^, 
and the continuum threshold fixed at sq = 2.3GeV^ as quoted in Ref. |Q. The sum rule is not 
optimized and 5^ = 0.3 is adopted. 

Here following Ref. p[, the sum rule is not optimized and the prediction for the nucleon 
tensor charge ^ 0.3 is obtained by simply taking = IGeV^. We see that the two sides 
of the sum rule does not match at all and they only coincide at the point = IGeV^, 
which is picked by hand. So, the prediction can be arbitrary as one may choose an arbitrary 
value for the Borel mass. 



IV. SUMMARY AND CONCLUSION 

In this paper, we have derived two new QCD sum rules for the nucleon tensor charge 
from the mixed correlator of spin- 1/2 and spin-3/2 interpolating fields. Both the usual 
formalism and the alternative formalism were adopted to handle the transitions between 
the ground state nucleon and excited states. We analyzed these new sum rules, along with 
the chiral-odd sum rule obtained from the usual 1/2 correlator, with respect to specific 
criteria for monitoring the contaminations arising from the transitions and continuum and 
the convergence of the OPE. An emphasis was put on revealing the validity and reliability 
of various sum rules. 

We found that in the usual formalism one of the two new sum rules is invalid and the 
other works only for certain values of x- The predictions of the latter for moderate x values 
(0.5 GeV2 <xa< 1-0 GeV^) are 
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(5m ~ 0.67- 1.06 , W~ 0.23 -0.24. (35) 

This shows that 6u has strong dependence on the vacuum tensor susceptibihty x and 6d is 
insensitive to x- The sum rule from the 1/2 correlator is valid for all the x values considered, 
and the predictions for the moderate x a-i's 

5m ~ 3.48 - 3.89 , M~ 0.12 -0.18. (36) 

It should be warned that the uncertainties quoted here only reflect the uncertainties from 
the X value. There are large uncertainties associated with the usual formalism arising 
from the treatment of the transitions, the spectral parameters in the mass sum rules, and 
the vacuum condensates and other susceptibilities. Without evaluating these uncertainties, 
we cannot make a critical comparison with the results from other models and approaches 
as these uncertainties may change these results significantly. Therefore, the above results 
should be treated only as benchmark (at most). Nevertheless, our results appear to support 
6u > 6d > 0. This behavior is different from that found in Rebs. p|-|TI| and the behavior of 
the quark spin structure of the nucleon (Am > 0, Ad < 0). 

The sum rules in the alternative formalism are invalid in most cases of interest. In 
particular, for isoscalar tensor charge all three sum rules are invalid. For isovector tensor 
charge, one of the two new sum rules is invalid for all the x values, and the other only works 
marginally for small x values and its predictions for the tensor charge are much smaller than 
those given in Eq. (|35|). The sum rule from the 1/2 correlator works only for large x values 
and its predictions are similar to those listed in Eq. (^). It is also found that the sum rule 
Eq. (^ with xa = -O.lSGeV^, studied in the previous works, is invalid. 

In conclusion, the QCD sum rules for the nucleon tensor charge work only in limited cases 
and their predictions are expected to have large uncertainties. It is important to examine 
the validity and reliability of the sum rules before the extraction can be made. Without such 
an examination, the extracted results are likely to be arbitrary and misleading. To improve 
the performances of the sum rules for the nucleon tensor charge, one should have more 
accurate evaluation of the vacuum susceptibilities and vacuum condensates, more precise 
knowledge of the mass sum-rule spectral parameters, higher order terms in the OPE, and 
better approaches for treating the transitions. 
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APPENDIX: 

In this appendix, we give some ingredients for calculating the QCD side as obtained from 
carrying out the OPE. Further details of the technique can be found in the literature. As 
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usual, we work to the leading order in perturbative theory and neglect the up and down 
current quark masses, the isospin breaking effect, and the strangeness contribution. It is 
convenient to calculate the Wilson coefficients of the OPE in coordinate space and then 
perform the Fourier transformation to momentum space. To first order in the external field, 
the quark propagator without background gluon field can be expressed as 



(T[gf(x)g?(0)]) 



X 



27648 



' 27r2 



* ^ij ^ab _ ill) 



5"' + ^ {qgs^ ■ Gq) S,, S'^' 



12 ''^ 192 ^' 



9q 



47r2 
9^ 1^ 



x^ 



288 



(gg) Z"'^ x" + 20 + 2 (fi: - ^, 



07 



(Al) 



where {«, j} and {a, 6} are Dirac and color indices, respectively, and the fifth term is the 
perturbative propagator in the presence of the external field. Note that the sign in front of 
^ differs from that given in Ref . p[ . Here we adopted the conventions of Itzykson and Zuber 
for e^fipx while Ref. |§] has not given the conventions exphcitly. 

In the fixed-point gauge, the quark propagator in the presence of background gluon field 
is given by 
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where the third and fourth terms arise from the presence of both the external and background 
gluon fields. The last three terms have been ignored in Ref. which gives rise to the 
differences in the last four terms in the sum rule Eq. (|27|). In the calculations of the terms 
proportional to ab, we have neglected the fourth and fifth terms of Eq. ([A^). We have found 
that the inclusion of such terms only make tiny refinement to the sum rules. 

We also need the following expression for the expectation value of the quark and gluon 
field product in the presence of the external field |TH 



ab 
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Now one can easily carry out the calculations of the correlators by expressing the corre- 
lators as iTBIl 
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with the quark propagator 5"^ (x, 0;^) = (T[g''(x)g (O)])^ given above. Since only linear 
response of the correlators is concerned, one only needs to keep the terms linear in the 
external field. 

Finally, we observe that the Borel transformation of ln(— g^)/g^ has not been carried out 
correctly in Ref. So, we give the correct formula 
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from which we get 
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It is then easy to show that the second term of Eq. (20) in Ref. |§] is not the correct Borel 
transform of the second term of Eq. (17). 



19 



REFERENCES 



R. L. Jaffe and X. Ji, Phys. Rev. Lett, 71, 2547 (1993). 

R. L. Jaffe, "S'pm, Twist and Hadron Structure in Deep Inelastic Processes" , Lectures 
presented at the International School of Nucleon Structure, Erice, Italy August 3-10, 
1995, Report No. MIT-CTP-2506, |hep-ph/ 96022^ , and references therein. 
J. Ralston and D. E. Soper, Nucl. Phys. B152, 109 (1979). 

RSC Collaboration, G. Bunce et al, Particle Word 3, 1 (1992); Proposal on Spin Physics 
using RHIC Polarized Collider, R5 (August 1992), approved October 1993. 
The HERMES proposal to HERA (1993). 

COMPASS Collaboration (G. Baum et al), Report No. CERN-SPSLC-96-14 (1996). 
H. X. He and X. Ji, Phys. Rev. D 52, 2960 (1995). 
H. X. He and X. Ji, Phys. Rev. D 54, 6897 (1996). 

H. -C. Kim, M. V. Polyakov, and K. Goeke, Phys. Rev. D 53, 4715 (1996); Phys. Lett. 
B387, 577 (1996). 

I. Schmidt and J. Sofferd, Report No. |hep-ph/97034Tl| . 



S. Aoki, M. Doui, T. Hatsuda, and Y. Kuramashi, Report No. hep-lat/9608115 . 

B. L. loffe, Nucl. Phys. B188, 317 (1981). 

D. B. Leinweber, Ann. Phys. (N.Y.) 254, 328 (1997). 

F. X. Lee, D. B. Leinweber, and X. Jin, Phys. Rev. D 55, 4066 (1997). 

B. L. loffe and A. V. Smilga, Nucl. Phys. B232, 109 (1984). 

X. Jin, Phys. Rev. D 55, 1693 (1997). 

Y. Chung, H. G. Dosch, M. Kremer, and D. Schall, Nucl. Phys. B197, 55 (1982). 

X. Jin and J. Tang, Report No. MIT-CTP-2601, |hep-ph/9701230| , Phys. Rev. D 56, 

(1997) in press. 

S. L. Wilson, J. Pasupathy, and C. B. Chiu, Phy. Rev. D 36, 1451 (1987). 
X. Jin and D. B. Leinweber, Phys. Rev. C 52, 3344 (1995); M. J. Iqbal, X. Jin, and D. B. 
Leinweber, Phys. Lett. B386, 55 (1996); R. J. Furnstahl, X. Jin, and D. B. Leinweber, 
Phys. Lett. B387, 253 (1996). 

V. M. Belyaev and A. Oganesian, Phys. Lett. B395, 307 (1997). 
B. L. loffe, Phys. At. Nucl. 58, 1408 (1995). 



20 



